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We obtain criteria for shock wave (SW) existence in relativistic 



O ■ magnetic hydrodynamics with with no suppositions about convex- 

ity of the equation of state. Method of derivation involves consid- 
eration of a continuous SW profile in presence of Landau-Lifshitz 



relativistic viscosity tensor with both non-zero viscosity coefficients 
r) and £. We point out that supposition of viscous profile existence 
with only one nonzero coefficient (rj=0) appears to be too restric- 
tive leading to losses of some physical solutions. 
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1 Introduction 

Relativistic shock waves (SW) arise in such powerful astrophysical phenom- 
ena as supernova explosions and gamma-ray bursts. Theoretical analysis of 
these processes involves criteria of existence and stability of discontinuous 
solutions that describe SW in superdense matter. Consideration of these cri- 
teria is complicated in case of general equation of state (EOS) (cf., e.g., the 
classical results 0] and the relativistic hydrodynamics [21 El)- 
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Note that in case of a normal fluid we deal with convex EOS (this means 
the convexity of Poisson adiabats); therefore the only condition is needed to 
study discontinuous solutions of hydrodynamical equations: the well known 
entropy growth criterion (see H] for classical hydrodynamics and jSJ E] in 
relativistic MHD). However, in case of a general EOS the convexity condition 
may be violated, and neither customary entropy criterion, nor the evolution- 
ary criterion jlj are not sufficient to single out physical solutions in a correct 
way. Moreover, the rarefaction shocks and the compression simple waves, as 
well as complicated configurations of shocks and simple waves moving in the 
same direction are possible. This situation is well-known in classical hydro- 
dynamics, it was first studied by H.Bethe |L. In the relativistic theory such 
anomalous equations of state arise, e.g., when dealing with a super-dense 
matter in the neighborhood of phase transitions (see, e.g., [21 IS])- 

One of the most effective methods to study the SW existence in case of 
the general EOS is investigation of the SW viscous profile. According to 
this method the generalized (discontinuous) solution is treated as a small 
viscosity limit of corresponding continuous solutions. The shock transition 
is admissible, if corresponding continuous solution (viscous profile) exists 
for any nonzero viscosity. In case of normal fluid (in the sense of Bethe and 
Weyl PP) the results of this method are the same as that of the evolutionarity 
criterion. In the relativistic hydrodynamics the conditions for viscous profile 
existence in case of the general EOS have been derived and studied by 
using the Landau-Liftshits viscosity term in relativistic energy-momentum 
tensor [4 . This term involves two viscosity coefficients £ and r\. 

In relativistic magnetic hydrodynamics (MHD) investigation of the 
SW viscous profile becomes more complicated. Therefore this problem has 
been first considered in a restricted version with one of the viscosity 
coefficients put equal to zero (rj = 0) under supposition that only one non- 
zero viscosity is sufficient to obtain the continuous SW profile. This was a 
technical supposition and it is not evident. At least, the results of 00 f° r 
rj = cannot be considered as necessary conditions. 

In the present paper we extend the results of to the case of arbi- 
trary ratio of positive viscosity coefficients and prove conditions for existence 
of the SW viscous profile under less restrictive requirements. We consider 
stationary viscous flows of relativistic fluid with infinite conductivity. These 
solutions describe the MHD shock structure, the existence of SW being con- 
sidered by means of corresponding continuous solutions with non-zero vis- 
cosity. Our treatment shows that we may relax the conditions of to 
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have a necessary and sufficient criteria. 

2 Basic equations 

The equations of motion of ideal relativistic fluid with infinite conductivity 
in magnetic field follow from the conservation laws involving the energy- 
momentum tensor jSJ EJ 

T» v = (p* + £*)u^u v - p*g^ u - ^-Wh v , (1) 

Air 

Vt? is the four velocity (Greek indexes run from to 3), the flat space-time 
metric g^ v = = diag(l, —1,-1,-1) is used for raising and lowering the 
indexes, = — |e Ma ^ 7 F a/ gM 7 is the magnetic field, e a/37<5 is the absolutely anti- 
symmetric symbol, F^ v is the tensor of electromagnetic field, p* = p + ^\h\ 2 , 
e* = e + ^r|/i| 2 , \h\ 2 = —h a h a > 0, /i is the magnetic permeability that is 
supposed to be constant; p is the pressure and e is the energy density (in 
the rest frame). We suppose EOS p = p(e,n) to be a sufficiently smooth 
function. 

Following the small viscosity method P3 El E] i n order to study the 
SW structure we introduce dissipation effects that smear out discontinuities. 
Similarly to [21 El HI E] > in case of relativistic problem we use the Landau- 
Liftshits viscosity tensor |lj 

V = V(u^ u + u u>ft - u^u a u^ a - u v u a u^ a ) + (f - 2r]/3)u a a (g f , u - u^u u ), 

the commas stand for derivatives. 

Now the fluid motion is constrained by equations of energy-momentum 
conservation 

baryon charge conservation 

fy(nu") = 0, (3) 
and one more equation follows from the Maxwell's equations El 

d^h u - u v W) = 0. (4) 
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The discontinuous solutions follow from these equations in the limit, when 
£ and r\ tend to zero. The questions is whether this limit depends on a relation 
between £ > and r\ > 0. 

The viscous profile of stationary SW may locally be represented in proper 
reference frame of the shock front by a stationary continuous solution depend- 
ing upon the only variable x; here t^ v — > and all the parameters of this 
viscous flow tend to constant values as x — > ±oo. 

Without loss of generality we suppose further that the limiting values 
of hydrodynamical parameters for x — >• — oo correspond to the state ahead 
of the shock (denoted further by index "0") and the values for x — > +oo 
correspond to the state behind the shock (denoted further by index "1"), 
then we have u 1 > behind and ahead of the shock. 

Because all the values in (JSJ),©,© depend only upon variable x, we have 
from these equations 

T^ + r^Tg, (5) 
u 1 ^ - h l u» = H» = u\ 0) hfa - h}^, (6) 

nu 1 = n(o)U( ). (7) 

As a result of t 1u — > for x — * ±oo, relations ©-(0) must be fulfilled 
for corresponding asymptotic values T^, n, obtained from continuous so- 
lutions of the system (J5j)- ()llj) . Similarly to classical hydrodynamics PQ we 
interpret the conditions for shock transition from the state u^, h^,n ,p 
(ahead of the shock) into the state u^,h^,ni,px (behind the shock), and 
as consequence of ©-© these states must satisfy 

ZW=T& (8) 

U (l)^(l) — ^\l) U \l) — H v = M(Q)/i( ) — fyo) u (0)> (9) 

n ( iju$ = 7i(o)U( )- (10) 
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3 Dynamical system for the shock structure 

In this section we use some of the results of jZHHJ- Suppose that equations 
P-fTTTj) are fulfilled. 

Definition. We say that shock transition u'Ly h^, Hq, po — > ufa, hf^s, nx,px 
has viscous profile if there is a continuous solution of having corre- 

sponding asymptotics for x — > — oo and x — > +00. 

We use the reference frame such that u 3 = and h 3 = 0; u 1 , /i 1 being 
normal components to the surface x = const, and u 2 , h 2 being the tangential 
components; «/L = 0. 

Following 0|H], due to (jHJ) we represent hP in terms of w 1 and u 2 : 

h? = —[H* -u»H a u a }. (11) 
Multiplying (JSJ) by and taking into account that r^ v u v = we have 

eV = 7$u M , (12) 

whence e and /i^ can be expressed through u l and u 2 . 

Taking into account the explicit form of r^ v we have from (j5J) for v = 1,2 



9m 1 



mil mil 

J (o) - J ' 



(13) 



1 



9m 2 
dx 



c ,V\ 1 
e+ 3j MM ^ 



rjil'2 rpYl 

1 (o) - 1 



(14) 



After eliminating of du 1 /dx from ()14j) with the help of (fT3*|) the second 
equation transforms to 



1 + 



du 2 
dx 



u l u 2 (f + ?7/3) 
f K) 2 (e + V3) 



mil mil 
1 (0) ~ J 



ml2 ml2 
J (0) ~ J 



It is convenient to introduce a new variable v — u Jl + (u l Y in (USD 
and (JHJ); this yields dynamical system with respect to u 1 and t> 



(15) 



where 
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Fi(y,v) = p- 



l + (u 



l\2 



+ 



(16) 



Fo(u\v) 



(T$u° - T^u 2 )u 2 u l - (fi/4ir)u 2 (H a u c 



(o) 



U 1[1 + ( U 1)2]! 



[1 + ( u i)a]3/a u i ^ J 

Continuous solutions of (|T5|) describe the SW structure. It is important 
to note that p(e,n) disappears from F 2 . 



4 Existence conditions of SW viscous profile. 

Let the state parameters u? s, h^,no, po ahead of the shock and ufo, h^,ni,pi 
behind the shock satisfy the conservation laws (|5|)- (fTU|) that relate hydrody- 
namic quantities on both sides of SW. We denote y = u l , yo = W( )> V 1 = m |tJ' 

Consider now the curves on (y, v) - plane where the right-hand sides 
(Unj), (fTTj) of system (fT3j) may change their signs. We denote Vi a locus of 
points (y,v) such that Fi(y,v) = and V 2 is a locus of points (y,v) such 
that F 2 (y,v) = 0. We suppose that Vi is represented by a smooth function 
v = Vi(y). From the results of |S] it follows that V 2 is represented by a 
single- valued function y = Y 2 (v); this function may be not monotonous so it 
is more convenient to use this function instead of the the inverse one. We 
suppose that "0" and "1" are connected by smooth components of Vi and 
V 2 (see, e.g., Figffl- FigH). 

We consider a part of (y, v) - plane between Vi and V 2 such that there 
are two intersection points "0" and "1" corresponding to the states ahead 
and behind the shock: v\ = Vi(yi), y x = Y 2 (vi); v = Vi(y ), y = Y 2 (v ), 
but the curves do not intersect between "0" and "1". 

The points (yo, Vi(yo)), (yi, Vi(yi)) are the rest points of system (|T5|). 

We shall consider the following conditions. 

A. The function v = V\(y) is a single- valued function on (yi,yo). 
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B. For all points of V 2 , v G (t>i,t>o) the following inequality is valid (cf. 

mm 

(y -y 1 )F 1 (Y 2 (v),v)<0. (18) 

Here we do not consider occurence of the Chapmen- Joguet points, where 
the left hand side of (jTSjl equals to zero but does not change its sign. This 
case may be studied by taking a corresponding limit in ([18)1 . 

C. We suppose that h x h 2 ^0 at the point "0". 

This is a technical requirement. Otherwise we deal with much more simple 
situation of parallel or perpendicular MHD SW; this case is not considered 
here. 

Let u s i, Uf, ua stand for the speeds of relativistic slow, fast and Alfven 
waves jniE]; v> s i, Uf being the roots of the polynomial Q(y), where 

Q(y) = (i - 4)(y 2 - u))tf - ul) = 

" (1 Cs)V V [° s+ A7r(p+e)) + 47r(p + e) [U} 

c l — (dp/de)s is the speed of sound. 

The relativistic Alfven speed ua is defined by the formula 

„. = (20) 

In [HE] one more parameter u* A has been introduced, which is the positive 
root of the equation R*(y) = 0, where 

R*(y) = (p+e)y 2 (l + y 2 )-£(h 1 ) 2 (21) 
The above velocities satisfy the inequalities jH] 

u sl < u A < u* A < Uf. (22) 

We introduce one more requirement (one can show that it is consistent 
with A,B). 

D. Either 

(DF): u 1 > Uf ahead of the shock at the point "0", 
or 
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(DS): u A > u 1 > u s i ahead of the shock at the point "0" . 

These inequalities correspond to the evolutionarity criteria of classical 
MHD for the velocities ahead of the shock [§]. The first inequality (DF) 
corresponds to the fast SW and the second one (DS) - to the slow SW. 

By means of (fTUj) - (fT2j) the variables p and e can be expressed in terms 
of the velocity components u l and u 2 . 

Lemma 1. If relations (J7j), - (J12j) are satisfied, then at the point 
"0" we have 



dS 



dS 



dv? 



0. 



where S is the entropy per baryon. 

The proof uses the thermodynamical relation TdS = pd(l/n) + d(e/n), 
where T is the temperature. The statement of the lemma is obtained after 
direct calculation of differentials dn and dp using ([TU ]) -([T2 |) . 

Using Lemma 1, after some calculations we have at the the point "0" 
(vo = 0): 



dF 1 



(p + e) 



Diu 1 ] 



(23) 



where 



D(y) 



c s V 



2cl + 



fi\h\' 



4tt (p + e) 



4vr (p + e) 



2\2l 



The other derivatives of the right hand sides (|lb|h (fT7j) of the dynamical 
system ([To} at the point "0" are 



dv 
8F 2 



du 1 
dF 2 



fi h l h 2 
h l h 2 



dv 



(vPYu 



(24) 
(25) 
(26) 



Taking into account these relations in the vicinity of the point "0" we 
have on the curve V 2 
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fito,V«»)) = ^^-«.). (27) 



On the curve V 2 we have 



47r«° J R*(w 1 )' ^ 

Lemma 2. In the case DF the rest point "1" of system (|15j) is a saddle 
point. In the case DS the rest point "0" is a saddle point. 
Proof. Direct calculation yields 



dFi dF 2 dF 1 dF 2 (p + e] 



2 



du 1 dv dv du l (u 1 ) 2 (u ^ 4 



In the case (DS) ahead of SW we have Q < at "0" (see flU) and 
this yields the required statement according to the properties of the saddle 
point [TU] . In the case (DF) analogous result at "1" can be checked directly; 
however it is easier to use the same relations as (|2"3 ]) -(|2lj |) in case of "0" 
by using the Lorents transformation that preserves ttL% and transforms the 
transversal velocity component u 2 ^ to zero at the point "1". 

Further the coefficients of Q(y) and R*(y) involving the magnetic field, 
energy density and baryon density are taken only at the point "0" ahead of 
the shock. 

Now we proceed to prove the existence of viscous profile. 

Consider first the case (DF) for y\ < y , and put for definiteness h l h 2 > Q 
at "0". Consideration of the opposite sign of h^h 2 is completely analogous. 

First of all we note that the condition (JTSJl guarantees that the trajectories 
of system (JT5j) in the plane (y = u 1 , v) cross the curve V2 from right to left 
(see Fig. QJ. 

On the curve v = Vi(y): 

, . . 47r (p + e) 2 Qlu 1 ) , . , , 

fi(v. *(»)) = - ; t! j;; 3 ) (y-y.)>° (29) 

in the vicinity of the point "0", y < yo. 

Correspondingly trajectories of system (|THj) cross the curve V\ bottom-up 
(Fig. [IJ. Evidently this is true not only for the neighborhood of "0", but for 
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all interval (yi,yo)', otherwise there must be additional rest points of system 
()15|) between "0" and "1" , which contradicts to our suppositions. 

Now we shall find out relative disposition of V\ and V%. Let at point "0" 



tg(a 1 ) 



dVi 



dy 



. . , dV 2 
tg(a 2 ) = — 



Using (f^ - fl^lj) we obtain that the ratio of tangents at "0" equals 



tg(a x ) _ 167r 2 (p + g )V) 2 _ 16vr 2 (p + e) 

tgj^j- W Q(n) + l- iflhW)2 R (« )D{u ) . (30) 

Then in case of fast SW {Q{yo) > 0) this ratio is > 1, that is \ x is above 

V 2 . 

Therefore, the phase curves only can leave the domain between Vi and 
V 2 (Fig. |TJ). Taking into account the sign of F 1; it is easy to see that inside 
this domain all the phase curves come out from the point "0" ; and there 
exists a phase curve of (jl5j) . that comes from "0" to "1" . Because "1" is 
a saddle point (Lemma 2), this phase curve is unique. This conclusion does 
not depend upon relation between (positive) viscosity coefficients £, r\. 

Now we proceed to the case (DS) (see Fig. |2J); it is now convenient to 
put h l h 2 < at "0" (the opposite case is completely analogous). The sign 
of (J2T?|) remains the same as in the case (DF) and so is the direction of the 
phase curves crossing Vi. The direction of the phase curves crossing V 2 also 
remains due to (fTSj) . Taking into account (DS) and we have Q(yo) < 
and according to (|3T?|) we see that V 2 is above V x in the neighborhood of 
"0". As distinct from the case (DF) here the phase curves can only enter 
the domain between "1" and "0". Similarly to previous consideration there 
is a unique phase curve of system (fTS|) . passing from "0" to "1"; this is a 
separatrix of saddle point " 0" . 

In the previous consideration we supposed that y\ < y ; this corresponds 
to usual compression SW. In case of anomalous EOS the rarefaction shocks 
are also possible P3E1IZIIH1- The condition (fTHj) is applicable in this case as 
well, and the consideration is completely similar. 

Therefore we proved that equations ©-(HUJ) have a unique continuous 
solution that connects the states "0" and "1". 

Theorem. Let the states u^, h^, n , p Q ahead of the shock and u^, 
h'ty rii, pi behind the shock satisfy the conservation equations (jHj)- (|10|) . If 
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the conditions (A-D) are satisfied, then the MHD shock transition " 0" — ^" 1" 
has a unique viscous profile satisfying equations (|5|)-(fTU]) . 

Note that the analogous criteria for existence of SW viscous profile ob- 
tained in [7JIB] appear to be too restrictive. These criteria have been obtained 
under condition that one of the viscosity coefficients equals to zero (77 = 0), 
and they rule out the shocks that satisfy the condition < u 1 < u* A at " 1" 
after the shock front. However these latter solutions are compatible with the 
criteria (A)-(D) of the present paper. The explanation of this inconsistency 
is as follows. If the function Y 2 (v) is monotonous and 77 — >- 0, then the phase 
curve of system (|T5|) that go from "0" to "1" tends to the curve V 2 . How- 
ever, it may happen that Y 2 (v) is not monotonous; this just corresponds to 
ua < u 1 < u* A at "1". In this case the above phase curve that begins at "0" 
snuggles down to V2 only on some segment, and the corresponding limiting 
solution for 77 —>■ has a discontinuity (see Fig. EJ). This explains why such 
solutions have been rejected in pHE], because the initial supposition of these 
papers was the existence of a regular viscous profile for 77 = . 

5 Discussion 

The SW existence conditions in relativistic MHD with general equation of 
state has been analyzed in |S] in case of 77 = in the Landau-Liftschits 
relativistic viscosity tensor. In the present paper we obtained the criteria 
for existence of viscous SW profile dealing with both nonzero viscosity coef- 
ficients (£ > 0, 77 > 0) in this viscosity tensor. If additional limitations on 
EOS (e.g., convexity) are absent, our criteria are more restrictive than, e.g., 
evolutionarity conditions jHj or any other conditions that involve character- 
istics of the fluid only at the initial and final states. This is evident because 
condition ()18|) must be valid for the whole interval between the states " 0" and 
"1" . This situation is analogous to ordinary (non-magnetic) hydrodynamics 
j2 El; in this case our criteria reduce to the criteria of these papers. 

On the other hand, the criteria of the present paper are less restrictive 
than that of 0|H] obtained in case of 77 = 0. This is because the supposition 
of existence of a regular viscous profile used in |H] does not always hold in 
case of 77 = (even if £ 7^ 0) and this rules out some physical solutions. It 
should be noted that such situation is specific just to relativistic MHD; this 
does not appear neither in non-relativistic case, nor in ordinary relativistic 
hydrodynamics. 
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Our criteria may be applied to arbitrary smooth EOS. However, we must 
note that the requirement for V\{y) to be a continuous (single- valued) func- 
tion is not trivial and may not be fulfilled in case of certain equations of state 
(cf., e.g., [1J). Though consideration of a viscous profile seems to be rather 
effective for investigation of SW existence and stability, this method may 
not work in case of complicated EOS (cf. remarks in ^T]) that require either 
modification of equations of the fluid motion or using additional physical 
information about solutions. 
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Figure 1: The phase trajectories of system (JT5j) in case of DF; EOS is p = 
e/3. The solid line (Vi) corresponds to v — V\{y), the dashed line (V 2 ) 
corresponds to y = Y 2 (v)\ the curve VP connecting "0" and "1" describes 
the viscous profile of the fast shock transition "0" — > "1". The arrows show 
the direction of the phase flow. 



13 



y 



Figure 2: The phase trajectories of system in case DS with the same EOS. 
In this case the curve Vj (solid) is below V 2 (dashed). The separatrix VP 
of the saddle point " 0" describing the the viscous SW profile of slow shock 
transition "0" -> "1" goes to the final state "1". 
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Figure 3: The example of viscous SW profiles for three ratios of in case 
of DF (fast SW) with non-monotonous dependence y = Y 2 {v); p = e/3. 
Disposition of curves Vi and V2 is as on Fig. ^ The curve "A" describes 
the profile with ratio rj/£ = 1, "B" corresponds to rj/£ = 0.1 and "C" - to 
i]/^ = 0.01. These curves corresponding to the smaller ratios snuggle down 
to V2 (after they go out from "0") on some segment and then jump to "1". 
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